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This article describes the performance of the all-digital data-transition tracking 
loop (DTTL) with coherent and noncoherent sampling using nonlinear theory. The 
effects of few samples per symbol and of non commensurate sampling and symbol 
rates are addressed and analyzed. Their impact on the probability density and 
variance of the phase error are quantified through computer simulations. 

It is shown that the performance of the all-digital DTTL approaches its analog 
counterpart when the sampling and symbol rates are noncommensurate (i.e., the 
number of samples per symbol is an irrational number). The loop signal-to-noise 
ratio (SNR) (inverse of phase error variance) degrades when the number of samples 
per symbol is an odd integer but degrades even further for even integers. 


I. Introduction 

In modern digital communication systems, analog-to- 
digital (A/D) conversion is performed as far toward the 
front end as possible using available technology. Usually, 
the received signal is amplified and then downconverted to 
the appropriate frequency for digital conversion. There- 
after, various system functions are performed digitally, in- 
cluding carrier, subcarrier, and symbol synchronization, as 
well as signal detection and decoding. Depending on the 
application, one can either sample the baseband signals 
(in-phase and quadrature) or sample the intermediate fre- 


quency (IF) signal. Furthermore, the sampling clock can 
be free running or controlled by the symbol synchroniza- 
tion loop. In the former case, the sampling rate is non- 
commensurate with the symbol rate. In the latter case, 
the sampling clock can be adjusted to obtain an integer 
number of samples per received symbol. All of these is- 
sues affect the final architecture and design of a receiver 
and influence the amount of cross-coupling among the var- 
ious loops. 

A receiver called the Advanced Receiver (ARX) has 
been developed for Deep Space Network (DSN) applica- 
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tions [1,2]. In the ARX, the signal is sampled at the 
IF, and the various tracking loops are implemented dig- 
itally, e.g., the “classical” analog integrate-and-dump fil- 
ters, which are typically part of the loop arms (in-phase 
and quadrature), are replaced by digital accumulators. 
This article investigates the performance of the all-digital 
data- transit ion tracking loop (DTTL), used for symbol 
synchronization, for any number of samples per symbol. 
In the previous version of the Advanced Receiver (ARX 
I) [1], the sampling was performed synchronously with the 
symbol rate, and an integer number of samples per sym- 
bol were available. In the current version of the Advanced 
Receiver (ARX II) [2], the sampling is performed asyn- 
chronously, and the sampling clock is fixed and indepen- 
dent of the symbol rate. At the highest desired data rate 
of 6.6 Msymbols/sec and with a fixed 20-MHz processing 
rate for the ARX II, only about 3 samples per symbol can 
be obtained. We are interested in the all-digital DTTL 
response and performance with a small number of samples 
per symbol. 

Some analytical results for the phase error variance of 
the analog DTTL were derived in [3], where the input is 
an analog signal and symbol and midphase detection were 
performed with analog integrate-and-dump filters. Later, 
the analysis was extended in [4], taking into account vari- 
ations of the equivalent noise spectrum with respect to 
normalized phase error. In this article, we extend the pre- 
vious results for the analog DTTL to the all-digital DTTL. 
We first note that there are two sampling models. One is 
to sample the signal instantaneously and the other is to ob- 
tain the sample by integration-and-dump (I & D) sampling 
of the signal. The instantaneous sampling technique can 
be used when the sampling rate is significantly higher than 
the symbol rate. The I & D sampling technique should 
be used when the number of samples per symbol is small 
[7]. In the absence of prefiltering and noise, the received 
symbol pulse is a perfect square wave, and instantaneous 
sampling provides voltage values of equal value. In the I 
& D case, all samples also have equal value except for the 
first one after a transition boundary of two symbols with 
different polarity. The all-digital DTTL can operate on 
either type of sample. 

To illustrate the differences between analog and all- 
digital DTTLs, we consider the noiseless case first. Note 
that if the input is an analog signal, the midphase inte- 
grator can produce a nonzero error voltage no matter how 
small the phase error is. Thus, a correction voltage can 
be generated at every symbol transition whenever a phase 
error exists. Therefore, the analog DTTL has infinite res- 
olution for phase detection. 


The all-digital DTTL, in contrast, has only finite res- 
olution for phase detection. This is illustrated in the fol- 
lowing: Suppose that there are an even number of samples 
per symbol. When a symbol transition occurs, the digi- 
tal midphase accumulator will produce a nonzero voltage 
only if the phase error causes sample slipping (assuming 
equal amplitude samples). As long as the phase error stays 
within a range of values that avoids sample slipping, the 
loop always generates a zero error signal. This range of 
undetectable phase errors accounts for the finite resolution 
of the all-digital DTTL. The more sampl es per symbol we 
use, the higher the resolution we can achieve and the closer 
the all-digital DTTL is to its analog counterpart. A key 
question is the impact of the all-digital DTTL’s finite res- 
olution on the S-curve and the phase error variance for few 
samples per symbol (say, three or four samples). 

Another issue in an all-digital implementation is the ef- 
fect of a noninteger number of samples per symbol. If the 
sampling clock is driven by the symbol synchronization 
loop, the number of samples per symbol can be made an 
exact even integer, which reduces the self-noise generated 
in the midphase accumulator, as will be discussed later. 
Under that sampling scenario, the sampling clock is con- 
stantly adapting as the data rate changes due to Doppler 
or other effects. One disadvantage of that scheme is that 
no fixed time base is available in the system. On the other 
hand, if the sampling clock is free running and is derived 
from a fixed frequency standard, the sampling period is 
fixed although the symbol rate may change. This may 
result in a noninteger number of samples per symbol. A 
model is derived in this article to analyze the performance 
of the all-digital DTTL for any sampling and symbol rates. 
In Section II, a nonlinear analysis of the loop is presented 
to handle all scenarios along with simulation results. The 
conclusion is presented in Section III. 


II. Analysis 

The performance of the all-digital data-transition track- 
ing loop with coherent and noncoherent sampling is ana- 
lyzed here. The block diagram of the all-digital DTTL is 
shown in Fig. 1. The input r(i) to the DTTL can be ob- 
tained by instantaneous sampling or by I k D sampling. 
In the subsequent derivation, closed form solutions will be 
obtained assuming samples of equal value, and results will 
be verified by simulation assuming samples of both equal 
and unequal value. 

Noncoherent sampling means that the sampling clock 
runs independently of the estimated symbol phase, i.e., 
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the sampling time is independent of the estimated symbol 
phase. This has minimal impact if there are many sam- 
ples per symbol. As the number of samples per symbol 
decreases, the S-curve becomes more and more coarse, and 
self-noise (to be explained later) increases. Our goal is to 
quantify the effects of a small number of samples per sym- 
bol on the phase error variance. A theory is presented for a 
first-order DTTL using Markov chains. The approach is to 
derive the S-curve and to solve the Chapman-Kolmogorov 
(C-K) equation to get the density function of the phase 
error. The closed-loop phase error variance and the degra- 
dation of the symbol detection can be evaluated from the 
phase error density function. 

To illustrate the phenomenon of self-noise, we consider 
a simple example shown in Fig. 2, where we have five sam- 
ples per symbol. We assume no thermal noise and perfect 
tracking at a particular moment. The output of the sym- 
bol transition detector is not zero because it either sums 
three samples from the first symbol and two samples from 
the second symbol [Fig. 2(a)] or sums two samples from 
the first symbol and three samples from the second sym- 
bol [Fig. 2(b)]. Notice that this situation occurs for ev- 
ery symbol interval as long as the loop maintains perfect 
tracking. The nonzero output of the loop filter will drag 
the loop away from the perfect tracking condition. 

In order to quantify this phenomenon, we first intro- 
duce three useful parameters. Let (3 denote the number 
of samples per symbol, which may not be an integer, and 
a(l) denote the offset of the first sample mark in a received 
symbol from the symbol boundary. By convention, a(&) is 
normalized and is measured as a fraction of the sampling 
interval. We observe that a(k) remains constant if (3 is an 
integer, and it varies from symbol to symbol if (3 is not an 
integer. Let us number the received symbols by 1, 2, 3, ... 
and denote the value of a(k) at the first symbol as a(l), 
which is referred to as the initial sampling offset. The 
values of a(fc) at the subsequent symbols, namely, o(2), 
c*(3), . . can be computed from (3 and a(l). The number 
of sample marks in a transition detection window and the 
number of sample marks in a symbol detection window 
are all functions of a(k). Thus, the output of the symbol 
detector and that of the transition detector fluctuate from 
symbol to symbol as a(k). This subject will be discussed 
later in more detail. 

Another important observation about the all-digital 
DTTL with a low number of samples per symbol is that it 
can have an irreducible phase error due to a finite number 
of samples per symbol. To illustrate this phenomenon, 
let us consider the example shown in Fig. 3, where ev- 


ery symbol contains four samples. We can see that as 
long as the estimated phase lies between 1 i and t 2i the 
error signal is always zero (or nearly zero if the received 
symbol does not have a perfect square waveform or if we 
have unequal amplitude samples) and the DTTL contin- 
ues tracking. However, we see that there still exists unre- 
solved phase ambiguity within the interval from t\ to t 2 . 
Mathematically, this phase ambiguity can be explained by 
a staircase S-curve. This phenomenon might have little 
effect on symbol detection performance if we use straight 
accumulation to detect the symbols. However, if we use 
weighted accumulation to detect the symbols, the phase 
ambiguity can introduce misweighting and thus degrade 
performance [5,6]. 

Before we proceed to the mathematical analysis, let us 
examine the all-digital DTTL block diagram again. The 
error signal accumulator between the loop filter F(z ) and 
the multiplier performs an averaging function so that the 
subsequent loop filter can operate at a slower speed. The 
loop bandwidth is determined primarily by the loop filter 
F(z). Thus, the presence of the accumulator is simply 
for hardware convenience. In the following analysis, we 
consider the DTTL without the error-signal accumulator. 

A. Mathematical Model 

Assuming that the carrier and subcarrier (if any) have 
been removed in an ideal fashion, the received baseband 
waveform is given by 

r(<) = %/S^ akp(t - kT) + n(t) (1) 

k 

where S is the signal power, T is the symbol time, n(t) is 
white Gaussian noise with one-sided power spectral den- 
sity No W/Hz, a* = ±1 represents the polarity of the 
jfcth symbol, and p(t) is the square-wave function having 
a value of 1 for 0 < i < T, and having a value of 0 else- 
where. The ith sample can be expressed after normalizing 
by 1 f\fS as 

r(i) = a* + n(i) (2) 

where we assume the sample is derived from the kth sym- 
bol, n(i) is a zero-mean Gaussian random variable with 
variance a 2 = N 0 /(2ST s ) y and T s is the sampling inter- 
val. Let the phase error A (in cycles) be defined as 
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where 9 is the actual received symbol phase and 9 is the 
estimated symbol phase. Note that A should have a value 
between -0.5 and 0.5. The error signal is affected by the 
locations of samples within their respective received sym- 
bols. In order to quantify this effect, we define a set of four 
A* integer-valued functions. They represent the number 
of sample marks contained in their respective intervals, as 
illustrated in Fig. 4. The output of the in-phase accu- 
mulator x(k) and the midphase accumulator y(k) can be 
expressed in terms of the A* functions, namely 


{ Ai (k)dk 4- Ao(fc)a.fc+i 4 n\(k) 4 n. 2 (A;) 
for A > 0 

Ai (k)ak 4 Ao(t)ujt-i 4- n\ (k) 4 n 2 (^) 
for A < 0 


x (k H- 1) 


r A\(k 4 1 ) a a- + 1 4- Ao {k 4 l)ajt+2 T n 3 (k) 
-f n 4 (A:) for A > 0 

A x(k 4 1 )a 1 T A 2 (& 4- 1 4* {k) 

< +n 4 (k) for A < 0 


y(k 4 1) = A3 (k)ak 4 A4(fc)a*+i 4 n 2 (k) 4- n 3 (k) (4) 


where rij(k),j = 1,2, 3,4 are zero-mean Gaussian random 
variables with variances 


Var{ni(fc)} 


Var{n 2 (fc)} 


r (Ai(fc) — A 3 (Jfc))<r 2 for A > 0 

< (Ai (*) + A 2 (jSr) + A 2 {k + 1) - A 3 (fc))<r 2 
, for A < 0 

f (A 3 (Ar) 4 A 2 (A:))^ 2 for A > 0 
\ (A 3 (k) - A 2 (k 4 l))^ 2 for A < 0 


Var{n 3 (&)} 


Var{n 4 (fc)} 


(A 4 (k) — A2 (k))<r 2 for A > 0 

(A 4 (k) 4 A 2 (k 4 1))<t 2 for A < 0 

(A 2 (i)4 A 1 (Jb4l) + A 2 (i4l) 
-A 4 (k))a 2 for A > 0 

(Ai(t 4 1) - A 4 (k))(T 2 for A < 0 


A !(*) 


[2/3 — a(k — 1)J — L(1 4 A)/? — a(k — 1)J 
for A > 0 

L(2 4 A)/? — a(k — 1)J — [p — at(k — 1)J 
for A < 0 


L(1 4 A)/? — ot{k)\ - [p-a(k)\ 
for A > 0 


A 2 {k) = i 

1 




IP - a(* - 

1)J-L(1 + A)/?- 

- a(k — 

1)J 


. for A < 0 



A 3 (*)= l/?-a(*)J- 

10 **-!)' 

- a(k ) j 


A 4 (i) = 

K 1+A+ 2- 

)/?-*(*)] - [/? 

- a(k) 

1 


with w denoting the width of the transition window and 
w < 1. As previously mentioned, a(k 4 1) is the sampling 
offset and is computed recursively from a(k) using 


a(* + 1) = VP - a(k) 4 1J -(p- a(k)) (7) 


and [t/J denotes the greatest integer strictly less than y 
(i.e., [4.2J — 4, [4J — 3). The derivation of the A * func- 
tions are similar and we illustrate only one example here, 
A 2 (i). To derive A 2 (fc), we use the beginning of the kth 
received symbol as the reference point. The number of 
samples in the fcth received symbol is [P — a(fc)J 4 1. The 
number of samples from the beginning of the fcth symbol to 
the end of the kth estimated symbol is [(I4A)/?— a(fc)J 41- 
Equation (5) follows by observing that the number of sam- 
ples from the end of the jfeth received symbol to the end 
of the kth estimated symbol is A 2 (&). The error signal 
e(k 4 1) is given by 


e(fc 4 1) — z{k 4 1 )y{k 4 1) (8) 


where z(ifc4l) denotes the output of the transition detector 
and sgn(z) denotes the “signum” function, i.e., 


„ , sgn [*(* + !)] - s S n i x ( k )] 

z(k + 1) = 5 


(9) 


(5) The conditional S-curve is defined by 


and the Aj functions are computed from 


g(\\a(k)) = E Btn {e{k+l)\\Mk)} ( 10 ) 
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where E s>n represents the conditional expectation with re- 
spect to both the signal (s) and the noise (n), assuming 
fixed A and a(k). Following mathematical manipulation 
similar to that in [4], we get 

^n{e(*+l)|A, <*(*)} = 

r 0.25[A 4 (*)(erf [r x {k)\ - erf [r 2 (fc)] + erf [r x (k + 1)] 

-|- erf [r 2 (fc + 1)]) - A 3 (fc)(erf [r x (k)] + erf [r 2 (ir)]) 

- r 3 (i)(A 3 (i) + A 2 (fc))(exp [— rj(fc)] + exp [— r|(fc)]) 
+ r 3 (i)(A 4 (fc) - A 2 (fc))(exp [— rj (k + 1)] 

+ exp [— r|(t + 1)])] for A > 0 
< 0.25[A 4 (ik)(erf [n(t + 1)] 4- erf [r 2 (i + 1)]) - A 3 (Ar) 
(erf [r x (k + 1)] - erf [r 2 (k + 1)] 

+ erf [r x (k)] + erf [r 2 (fc)]) 

- r 3 (fc)(A 3 (*) - A 2 (^))(exp [-rj(4)] + exp [— rf(*)]) 
4- r 3 (fc)(A 4 (fc) + A 2 (fc))(exp [-r\ (k + 1)] 

k exp [— r?(fc + 1)])] for A < 0 

(11) 


We observe that £f(A|a(fc)) is the (unconditional) S-curve 
if a(ife) is a constant. If a(k) changes from symbol to sym- 
bol, its effect will be smoothed in the loop as long as the 
loop time constant (inverse of loop bandwidth) is larger 
than MT. M denotes the number of distinct a(k ) y s and 
is discussed in the subsequent section. Therefore, the S- 
curve is obtained by averaging Eq. (11) over the M values 
of a(Jfc), which is determined by the initial sampling offset 
and the number of samples per symbol /?, that is, 

1 M 

ff(A) = tt5>(A|«(*)) (13) 

1 Jfc=l 

The S-curve #(A) is a function of A, /?, and a(k). In 
general, j3 can be any real number. When (3 is an irrational 
number, then a(k ) is nonperiodic. However, when j3 is a 
rational number, then a(k) becomes periodic and assumes 
a finite number of possible values. For example, if (3 is of 
the form X.0000 . . where X is any nonzero integer, then 
a(k) will have only one value all the time; if (3 is of the 
form X.X000 . . ., then a(fc) will have at most 10 different 
values; if j3 is of the form X. XX 000 . . ., then a(k) will have 
at most 100 different values, and so forth. For example, 
consider j3 — 4.1 and suppose the initial sampling offset 
a(l) = 0.7. Clearly, a (2) = 0.6, a(3) = 0.5, a(4) = 0.4, 
a(5) = 0.3, a(6) = 0.2, a( 7) = 0.1, a( 8) = 0, a(9) = 0.9, 
a(10) = 0.8, o:( 11) = 0.7 .... The more values c*(£) takes, 
the smoother will be the S-curve of the tracking loop. 


where 


R s (A x (k) - A 2 {k)Y 
(3(A x (k) + A 2 (i)) 


For a large /?, ri(k + 1) = r,*(fc) and i = 1,2,3; for 
A > 0, Ai(fc) -+ (1 - A)/?, A 2 (k) -* A/?, A 3 (k ) (w/2 

- A)/?, A 4 (Jfe) — ♦ ( w/2 + A)/?; and, as a result, ri(k) —* (1 

— 2X)y/R 7, r 2 (k) — ♦ \/W s , r 3 (k) — ► 1/ y/n W s . Substituting 
back in Eq. (11) and simplifying, one obtains 



g..n{^e (A)} =A erf( N /fl7(i- 2 A)) 

- l[u; - 2A] [erf (x/^T) 

-erf( N /fl7(l-2A))] (14) 

which agrees with the S-curve derived by Simon [4] for the 
analog loop. For A < 0, similar steps can be taken to show 
that the resulting S-curve also agrees with the analog case, 
with slope at the origin given by 


208 


m = P ^erf {\fR,) - ( 15 ) 

Figure 5 depicts the theoretical and simulated S-curves as 
a function of the phase error for various values of /?, namely 
4, 5, 10, 4.5, and 4.74, with w = 1. Note the staircase 
shape of the S-curve for all /?, with step size proportional 
to 1 /M/3. In this case, the phase detector is insensitive to 
variations of the input phase which occur in the flat regions 
of the S-curve. For small integers, the steps are large, while 
for large integers, the steps are small. In general, the S- 
curve is not symmetric and has a bias (i.e., #(0) ^ 0) 
proportional to 1/M/?, except for the case when /? is an 
exact even integer. The sign of the asymmetry depends on 
a(l), which was set to 0.5 in Fig. 5(e). As /? approaches 
an irrational number, the S-curve becomes smooth due to 
the averaging over a large set of values and it converges 
to the S-curve encountered in analog systems. In all these 
figures, the S-curve has a zero slope at the origin, which 
prevents a linear analysis because it requires a nonzero 
slope. Also depicted in Fig. 5(a) are two cases: when the 
signal and noise (s + n) are filtered prior to the loop and 
when only the signal is filtered but not the noise. The first 
scenario represents the case where the main contribution 
of the filtering is occurring at the receiver, while the latter 
scenario is representative of filtering at the transmitter. 
The filter used in the simulations is a first-order lowpass 
with transfer function H(z) = ho/ (l + hiz” 1 ). Note from 
Fig. 5(a) that the filtering introduces a bias and a slight 
asymmetry in the S-curve which is mainly due to the signal 
filtering. The asymmetry in the S-curve can be reduced 
by using a raised cosine or square-root raised cosine filter 
[7] that results in a symmetrical pulse shape around the 
data transition point. In the remaining figures, simulations 
with both signal and noise filtering are shown to assess the 
impact of imperfect square pulses. 

One way to linearize the loop and, hence, to smooth the 
shape of the S-curve is to randomly shift by an amount the 
position of the window relative to the symbol transition 
point. The effect of this random back and forth shifting 
is to produce on the average the same amount of samples 
in both halves of the symbol transition detector, so that 
/? would appear to have an irrational value to the mid- 
phase detector. The performance of the DTTL with such 
a detector can be the topic of future work. 

We will proceed to perform a nonlinear analysis of the 
loop using the C-K equation. However, we will need the 
open-loop variance of the error signal at the input to the 
loop filter to compute the transitional probability density. 
The second moment of the error signal is given by 


e2(A) = ^, n |^f;e(^A) 

1 M „ Af-1 

= T7T / E s , n {e(k,Xy} + w - -J2 E.,n MM) 

1 k = 1 1 k - 1 

x e(k + 1, A)} 

1 M 0 Af-1 

= m S h (t) + + B(m - i) S 

x [t> 3 (i) + v ^(k) + t>5(fc)] (16) 

where the variables Vi(k ), i = 1,2, 3,4, 5 are defined in 
the Appendix. The open-loop variance is then computed 
using 

<7 e 2 (A) = e 2 (A)- fl 2 (A) (17) 

The open loop variance a^( A) is depicted in Fig. 6 for 
various values of /? as a function of A. Again, note the 
staircase shape of the variance which approaches a smooth 
function for /? — 4.74. The analysis agrees very well with 
the simulation points for all values of /?. 

B. Probability Density Function of the Phase Error 

When the number of samples per symbol, /?, is low and 
the decimal part of /? has only a few nonzero digits, the 
S-curve takes a coarse staircase shape, as discussed previ- 
ously. In this case, the loop behaves nonlinearly, and the 
phase jitter of the loop cannot be predicted from the loop’s 
noise-equivalent bandwidth. However, when the statistics 
of the driving noise process are known, the C-K equa- 
tion permits one to derive the probability density function 
(pdf) of the closed-loop phase jitter even for a highly non- 
linear loop. Once this pdf is found, all the moments of 
the phase jitter process can be computed, and the noise 
performance of the loop can be predicted. 

In general, an all-digital phase-locked loop can be de- 
scribed by a stochastic difference equation of the form 

A (k) = u(k) - [g(X(k)) + n(k)]N(z)F(z) (18) 

where A (k) is the phase error at the instant k } u(k) is 
the input signal sequence, F(z) and N(z) are the dis- 
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Crete transfer functions of the loop filter and the numer- 
ically controlled oscillator (NCO), respectively, the func- 
tion < 7 (A(&)) represents the nonlinearity of the phase de- 
tector, and n(k) is the open-loop noise process. In digi- 
tal phase-locked loops, this noise process consists of both 
thermal and quantization noise. For the all-digital DTTL 
under consideration, the noise n(fc) is zero-mean with 
the variance given by Eq. (17). In general, the product 
F(z)N(z) is of the form 


a = [ai a 2 • • • ajv]* , and b = [iq b 2 ■ • * b^Y (the superscript 
t denotes transpose). The need to rewrite the original 
stochastic difference in vector form stems from the fact 
that the original equation does not represent a Markov 
process while its new form is a vector Markov process. 
For a vector Markov process, the C-K equation relates the 
probability density function of the state vector x at time 
(jfe 4- 1), pk+i(x(k + 1)), to its probability density function 
at time k through the following integral: 


F(z)N(z) 


bo 4* b\z 1 + ...-b 6j\rz N 
do 4 - a\z~ l + . . . + aj^z~ N 


( 19 ) 


Here in the context of the stochastic difference equations, 
z~ l should be thought of as the unit delay operator. Using 
Eq. (19), Eq. (18) can be rewritten as 


A (k) — [a Q u(k) + b 0 {g{X(k)) + n{k))] 

4- z~ l [aiu(k) - a\X(k) - big(X(k)) - bin(k)] 


P*+i(x(* + 1) = *l*o )=/•••/ q\x(k + 1) = x\x(k) = 

J — oo j 


~oo . 
N—dim 


v]Pk{x(k) = v\ Xjo) dv 1 - * ■ dv N 

(23) 


where z# is the initial condition vector and q(-) is the tran- 
sition probability density function. 


-f z 2 [a 2 u(k) - a 2 X (k) - b 2 g(X(k)) - b 2 n(k)] 


+ 


-f z N [a N u(k) - a^X(k) - b^g(X(k)) - b N n(k)] 

( 20 ) 

Defining x\ as the row with z -1 , x 2 as the row with z“ 2 , 
etc., the above equation can now be expressed in terms of 
the following state and output equations: 

x(k 4* 1) = A x(k) — bg(X(k)) + au(k) - bn(k) 


We will now focus on solving the above equation for the 
steady-state case (large £) and a periodic phase detector. 
In steady state, the initial condition vector Xq is washed 
out and can be dropped from the C-K equation, p* + i( ) — ► 

Pk(') = p( ) j so that the index k can be dropped also. With 
a periodic phase detector <;(*), each state variable takes 
values only between 0.5 < Xi < 0.5 (i = 1 ...N). This 
requires the “folding” (or collapsing) of the transitional pdf 
q[-] to the bounded region. We define this new transitional 
pdf as £[■] and rewrite the C-K as follows: 


p{y = y)= f ••• /<;[y(fc + i) = 2/|y(k) 

J — 0.5 J 


N — dim 


\(k) = xi(k) + [a 0 u(k) 4- b Q g(X(k)) 4- b 0 n(k)] 


— 1 ^] P( y — w)dwi ■ ■ ■ dujj\ ( 24 ) 


( 21 ) 


where x_ is the state vector, [x\ x 2 ■ ■ • 
sition matrix 

’ — ai 10 0 

— a 2 0 10- 

• 0 0 ■ ■ 

. — a j v 0 0 0 0_ 


A is the tran- 


( 22 ) 


where the primed p represents the new pdf restricted to 
the {—0.5, 0.5} region. Even in this simplified form of the 
C-K equation, closed form solutions can be obtained only 
for a few special cases. In general, the C-K equation has 
to be solved numerically, and solutions for N larger than 
1 can become very computationally demanding. 

In order to proceed, we will restrict our attention to the 
case of N = 1 with F(z) — bi and N(z) — z ~ l / 1 — z _1 
(ideal summer), which corresponds to a type-I loop with no 


210 


computational delay. Rewriting Eq. (19) for this restricted 
case, we get 


where a] = b\c r*(A), <r 2 e ( A) is the variance given by 
Eq. (17) at A - Xj t and fij = -a^Xj - b l g(x j ). 


F(z)N(z) 


b\z 1 
a 0 + a x z~ l 


(25) 


With N = 1, the variables in Eq. (20) are ao = 1 and a\ = 
— 1 with all others being zero. Since we are not considering 
any static phase error, we set the input signal Uk — ► 0 so 
that the loop is driven by noise only. Next, we discretize 
the continuous variable y into L + 1 values and, thus, 
approximate the Markov process by a Markov chain [8]. 
The larger the value of L , the better is the approximation. 
The continuous pdf on both sides of the C-K equation can 
be replaced by discrete probabilities P(yi) — 5p(y = yi) 
where 6 = 1/L, and the transition pdf £[i,j] is replaced by 
transition probability mass distribution Q[iJ] = $ 2 £[i, j]. 
With this substitution, the C-K equation becomes 


P(y — yi) — ^ Q[y(k + 1) — Vi\ vW — w j] P{y — w j) 

(26) 


i = " 31 


However, the Gaussian assumption becomes less accu- 
rate when the quantization noise starts to dominate over 
the thermal noise (the high-SNR case). Our simulations 
have shown that, at high SNRs, the pdf on the phase er- 
ror process becomes highly irregular and difficult to de- 
scribe mathematically. So a trapezoidal or a uniform dis- 
tribution function with mean and variance as used for the 
Gaussian gives results close to the ones obtained by simu- 
lation. In our computations, the singular value decompo 
sition (SVD) algorithm was used, and L was set to 1100 
to achieve acceptable resolution at a high SNR. The com- 
puted and simulated pdfs are depicted in Fig. 7 for the 
case of 0 = 4,5,10,4.5, and 4.74. When 0 = 4, the pdf 
exhibits a significant flat region, as expected. However, 
for 0 = 5, the flat region disappears due to the effect of 
the odd number of samples per symbol. The simulations 
were carried out with one-sided loop bandwidth Bi set to 
0.01 Hz and a symbol rate of 1 symbol/sec. For 0 = 10, 
the flat region is present again but is significantly reduced 
compared to the case of 0 = 4. For 0 = 4.5 or 4.74, the 
pdf is a smooth function, as one would expect. Note the 
nonzero mean in Fig. 7(e), which is due to the asymmetric 
S-curve in Fig. 5(e). 


where Xi — i6 } i = — i/2, • • ■ L/2. The above equation has 
to be true for all i. Let £(■) be the (L T 1) dimensional 
vector with elements P(y = yi) and Q be the (L + 1) x 
(L+l) matrix with elements Q[i, j]\ then the C-K equation 
can be written in the compact form P_ = QP, which can 
be solved by various techniques used in systems of linear 
equations [9]. For example, one can write the above matrix 
equation as 


[Q - I]P = 0 (27) 


with I being the identity matrix and 0 the zero vector. 
This matrix Q — I will have at most (L + 1) distinct eigen- 
values, and the desired solution, the P vector, will be the 
eigenvector corresponding to the zero eigenvalue. The last 
step in solving for P is to assume a suitable transition 
probability density function which will generate the ele- 
ments Q[iJ] of the matrix Q. When thermal noise domi- 
nates, the noise process n(Jc) can be assumed to be Gaus- 
sian, for which 


Q[iJ] = 




ex p 


(*; + /- //y) 2 


2 ( 7 ? 


(28) 


The corresponding variances for all cases are shown in 
Table 1. In Table 1(a), the loop SNR (inverse of tracking 
variance) is higher at 0 = 5 than at 0 = 4 and 10 due 
to the averaging process over the various offsets, a(k)' s. 
Note that the simulations and predictions agree well in all 
cases, with the largest deviation, of about 0.7 dB, obtained 
with 0 = 5. Table 1(b) depicts the effect of the self-noise 
which is dominant at high symbol SNRs (such as 25 dB). 
The model predicts the contribution of the self-noise very 
accurately for integer as well as noninteger 0. This noise 
is the limiting factor in the tracking performance in any 
all-digital DTTL implementation and is nonexistent in an 
analog design. 


III. Conclusion 

Nonlinear analysis of the all-digital DTTL is used to 
quantify its performance as a function of the number of 
samples per symbol 0. It is shown that the probability 
density function of the closed-loop phase error can vary 
significantly depending upon the number of samples per 
symbol, 0, and the symbol SNR. 

The performance of the all-digital DTTL approaches its 
analog counterpart as 0 increases and the sampling and 
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symbol rates are noncommensurate (i.e., ft is irrational). 
The loop SNR (inverse of phase error variance) degrades 
when P is an odd integer and degrades even further when 


ft is an even integer. In general, the S-curve has a bias 
proportional to 1/(M/?), but for an even-integer ft the bias 
goes to zero. 
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Table 1 . Predicted and simulated variances: (a) white noise 
dominated region and (b) self-noise dominated region. 


0 

Symbol 
SNR, dB 

Variance, 

simulation 

Variance, 

theory 

4.0 

3 

0.00635 

0.006324 

5.0 

3 

0.000954 

0.00114 

10.0 

3 

0.00234 

0.002434 

10.0 

10 

0.000925 

0.000972 

4.5 

3 

0.00235 

0.0026 

4.74 

3 

0.00253 

0.002720 


(b) 


P 

Symbol 
SNR, dB 

Variance, 

simulation 

Variance, 

theory 

4.0 

25 

0.004448 

0.004578 

5.0 

25 

0.000028 

0.000030 

10.0 

25 

0.000910 

0.000972 

4.5 

25 

0.000021 

0.000024 

4.75 

25 

0.002633 

0.002750 



Fig. 1. Block diagram of the all-digital DTTL. 



I 


Fig. 2. The effect of an odd number of samples per symbol: 
(a) three samples from the first symbol and two samples 
from the second symbol and (b) two samples from the first 
symbol and three samples from the second symbol. 
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Fig. 3. The phase ambiguity phenomenon. 
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Fig. 4. Definitions of A/ functions: (a) A > 0 and (b) A < 0. 
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Fig. 5. S curves for various /3’s: (a) (3 = 4, symbol SNR = 3 dB, (b) ( 3 = 5 , symbol SNR = 3 dB, (c) (3 
symbol SNR = 3 dB, (d) (3 = 10, symbol SNR = 10 dB, (e) (3 = 4.5, symbol SNR = 3 dB, and (f) /3 = 
symbol SNR = 3 dB. 











0.7 




















Fig. 7. Probability density function p( A): (a) ( 3 = 4 , symbol SNR = 3 dB, (b) /3 — 5, symbol SNR = 3 dB, 
(c) 0 = 10, symbol SNR = 3 dB, (d) (3 = 10, symbol SNR = 10 dB, (e) 0 = 4.5, symbol SNR = 3 dB, and 
(f) 0 = 4.74, symbol SNR = 3 dB. 



















Appendix 

Definition of Variables v,{k), v 2 (k), v 3 {k), v 4 {k), and v 5 {k) 


vi{k) = 


A 2 (fc) + A 2 (fc) + 


m aW + MM 

2R, 


x [4 + (erf [ri(fc)] - erf [r 2 (fc)])(erf [r^fc + 1)] 4- erf [r 2 (k + 1)])]^ 

- 2A 3 (&)A 4 (£)(erf [rj (*)] + erf [r 2 (fc)])(erf [ri(k + 1)] + erf [r 2 (fc + 1)]) 
+ (A 3 (fc) + A 2 (fc)) 2 r 3 (fc) ^(erf [ri(k + 1)] + erf [r 2 (fc + 1)]) 

x (- Mtl + tw ” 1 * [_rf(t + 1)1 + " p [_r5( * + 1)] )) 

- (A 4 (fc) - A 2 (k)) 2 r 3 (k + 1) ^(erf [ri(i)] - erf [r 2 (fc)]) 


x 


Ai(t + 1)- A 2 (fc + f) 
A + 1) + A 2 (k + 1) 


(exp [-r 2 (fc + 1)] + exp [~rl(k + 



for A > 0 


(A-l) 


vi(k) = 


&i(k) + A \{k) + 


/?[A 3 (Jfc) + A 4 (*)] 


2R, 


x [4 + (erf [ri(Ar)] - erf [r 2 (fc)])(erf [^(fc + 1)] + erf [r 2 (fc + 1)])]^ 


- 2A 3 (fc)A 4 (fc)(erf [r^fc)] + erf [r 2 (fc)])(erf [ri(fc + 1)] + erf [r 2 (k + 1)]) 
+ (A 3 (fc) - A 2 (fc)) 2 r 3 (fc) ^(erf [ri(k + 1)] + erf [r 2 (fc + 1)]) 

*(- £$ i" r i (l+ 1) i +exp [—•!(*+ i>i)) 

- (A 4 (ife) + A 2 (£)) 2 r 3 (fc + 1) ^(erf [ri(k)] - erf [r 2 (A:)]) 


X 


/ Ai [k + 1) — A2(fc + 1) 
\ Ai(fc -j- 1) + A2(^ + 1) 


exp [~rl(k + 1)] + exp [-rf 



for A < 0 


(A-2) 
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v 2 (k) = 2A 3 (k)(r 3 (k)[A 3 (k) + A 2 (fc)](erf [n(fc + 1)] + erf [r 2 (fc + 1)]) 

x (exp [— rl(k + 1)] — exp [— r\{k + 1)]) — r 3 (k + l)[A 4 (fc) — A 2 (fc)] 
x (erf [ri(k)} + erf [r 2 (fc)])(exp [-r 2 (k + 1)] + exp [-r\{k + 1)])) 

+ 2A 4 (k)(-r 3 {k)[A 3 {k) + A 2 (fc)](erf [r^k + 1)] + erf [r 2 (k + 1)]) 
x (exp [-r? (k + 1)] + exp [-rj(k + 1)]) + r 3 (k + l)[A 4 (fc) - A 2 (fc)] 

x (erf [n(k)] - erf [r 2 (fc)])(exp [-rj(fc + 1)] + exp [-r^fc + 1)])) for A > 0 (A-3) 


v 2 (k) = 2A 3 (Jfe)(r 3 (Jfc)[A 3 (Jfc) + A 2 (fc)](erf [ ri (k + 1)] - erf [r 2 (k + 1)]) 

x (exp [— r\(k + 1)] + exp [— r 2 {k + 1)]) — r 3 (k + l)[A 4 (Ar) — A 2 (&)] 
x (erf [rj(£)] + erf [r 2 (fc)])(exp [-ri(fc + 1)] + exp [-^(fc + 1)])) 

+ 2A 4 (fc)(-r 3 (fc)[A 3 (fc) + A 2 (fc)](erf [ri(k + 1)] + erf [r 2 (k + 1)]). 
x (exp [~r\(k + 1)] + exp [~r 2 (k + 1)]) + r 3 (k + l)[A 4 (fc) - A 2 (fc)] 

x (erf [ri(ife)] + erf [r 2 (fc)])(exp [-r^fc + 1)] - exp [-r 2 2 {k + 1)])) for A < 0 (A-4) 


v 3 (k) = A 3 (k)A 3 (k + l)[(erf [n(fc)] + erf [r 2 (£)])(erf [n(k + 1)] + erf [r 2 (k + 1)])] 

- (A 3 (fc)A 4 (fc + l)(erf [n(k)] + erf [r 2 (fc)])((erf [n(k + 1)] - erf [r 2 (k + 1)]) 

+ (erf + 2)] + erf [r 2 (k + 2)]))) 

- A 3 (k + l)A 4 (fc)(4 + [(erf [ri(k)] - erf [r 2 (fc)])(erf [r^k + 1)] 

+ erf [r 2 (k + 1)])] + (erf [n(k + 1)] - erf [r 2 (k + 1)]) 

x (erf [ri(k + 2)] + erf [r 2 (k + 2)])) 

+ A 4 (fc)A 4 (fc + l)((erf [r^k)] - erf [r 2 (*)])((erf [r^k + 1)] - erf [r 2 {k + 1)]) 

+ (erf [ri(k + 2)] + erf [r 2 (k + 2)])) + (erf [r^k + 1)] + erf [r 2 (k + 1)]) 
x (erf [n(* + 2)] + erf [r 2 (* + 2)])) for A > 0 (A-5) 
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v 3 (k) = A 4 (k)A 4 (k + 1 ) [(erf [ri(k + 1)] + erf [r 2 (fc + 1)]) 
x (erf [rj(fc + 2)] + erf [r 2 (k + 2)])] 

- (A 3 (k)A 4 (k + l)(erf [r t (fc + 2)] + erf [r 2 (k + 2)])((erf [ri(fc)] + erf [r 2 (k)}) 
+ (erf [n(*+l)]-erf [r 2 (fc + 1)]))) 

- A 3 (k + l)A 4 (fc)(4 4- [(erf [ri(fc)] + erf [r 2 (A;)])(erf [ri(k + 1)] 

- erf [r 2 (k + 1)])] + (erf [r x (k + 2)] - erf [r 2 (fc + 2)]) 

x (erf [ri(fc + 1)] + erf [r 2 (fc + 1)])) 

+ A 3 (k)A 3 (k + l)((erf [»*i(*)] + erf [r 2 (fc)])((erf [r 4 (k + 2)] - erf [r 2 (fc + 2)]) 
+ (erf [ri(k + 1)] + erf [r 2 (k + 1)])) + (erf [ri(k + 1)] - erf [r 2 (k + 1)]) 
x (erf [ri(fc + 2)] — erf [r 2 (fc + 2)])) for A < 0 


v 4 (k) = (r 3 (k)[A 3 (k) + A 2 {k)]r 3 (k + l)[A 3 (fc + 1) + A 2 {k + l)](exp [~r^(k + 1)] 

+ exp [— r 2 (fc + l)])(exp [-r\{k + 1)] + exp [-r\(k + 1)])) 

+ (A 3 (k)[A 4 (k) - A 2 (k)]r 3 (k + l)[A 4 (fc + 1) - A 2 (k + 1)] 
x (exp [~rl(k + 1)] + exp [~rl(k + l)])(exp [-r?( k + 2)] + exp [-r|(fc + 2)])) 
- (r 3 (Ar)[A 3 (Jfc) + A 2 (*)]r 3 (* + 1)[A 4 (Jb + 1) - A 2 (* + l)](exp [-r?(i + 1)] 

+ exp [~r%(k + l)])(exp [-r?(fc + 2)] + exp [~rl(k + 2)])) 


- (erf [ ri(k )] - erf [r 2 (*)])t 


/ /?(A 4 (*) - A 2 (fc)) 4 

ttR, 


(Ai(k + 1) — A 2 (fc + 1)) 


(A 4 (Jfc) - A 2 (Jfc) + A 3 (k + 1) - A 2 (k + l)) 1 ' 


■ exp 


-(Ai(A: + 1) - A 2 (k + l)) 2 -^ 


[A 3 (k) + A 2 (*)] 


1 + 


^ 3 ^ + 1 ) + + 1 ) 

lA 4 (t) - Aa(*)J 


+ 


( 


(A t (ib + 1) + A a (t + 1)) 


(A 4 (k) - A 2 (k) + A 3 (Jb + 1) - A 2 (k + l)) 1 ; 


■ exp 


[A 4 (k) - A 2 (fc)] 


-(A^ + O + A^ + l)) 2 ^ 


(A-6) 


221 



[A 3 (*) + Aa(*)] 


1 + 


A 3 (£ 4 1 ) 4 A 2 (fc + 1 )] ] [A^fc) — A 2 (£)] 

;a 4 (*)- a 2 (*)] 


+ (erf [ri(k + 2 )] -f erf [r 2 (k -f 2 )])* 


! P(A 4 (k) - A 2 (k)Y 
xR s 


x - 


(Ai(fc 4- 1) — A 2 (k + 1 )) 


(A 4 (Ar) — A 2 (k) 4- A 3 (& + 1) — A 2 (k + l )) 1 - 5 


exp 


— (Ai(fc + 1) - A 2 (fc 4 1 )) 


2 ^ 

~T 




1 

1 

n 

[Aa(*) + A 2 (*)] 

1 + J 

A 3 (fc 4* 1) 4~ A 2 {k 4- 1) 

[A 4 (Jb) — A 2 (*)] 


IA 4 (*) - A 2 (jb)] J 




+ 


(Ai(& 4- 1) 4* A 2 {k 4~ 1)) 


(A 4 (fc) - A 2 (Jfc) 4 A 3 (* + 1) - A 2 (* 4 l)) 1 5 
1 


exp 


-(Aitfc+O + AaOfc+l)) 2 ^ 


[A 1 A (Mlfl 1 [ a 3(^ + 1) + A 2 (fc + 1)] 1 [A 4 (fc) - A 2 (fc)] 

[Aa( ^ + A2 ^ )J i 1 + [A 4 (t)-A 2 (ibjJ J 


(A-7) 


v 5 (k) = r 3 (fc) [A 3 (k) + A 2 (Ar)](cxp [-r?(fc + 1)] + exp [-r|(fc + 1)])(A 3 (£ + l)(erf [r^k + 1 )] + erf [r 2 (fc + 1)]) 

- A 4 (k + l)[(erf [r x (fc + 1)] - erf [r 2 (k + 1)]) + (erf [r x (k -f 2)] + erf [r 2 (fc + 2)])]) + A 3 (fc + 1) 

x [A 3 (fc + L) + A 2 (fc + l)](A 3 (i')(erf [r x (k)] + erf [r 2 (fc)])(exp [-r x (A: + 1)] + exp [-r? 2 (k + 1)]) 

- A 4 (fc)[(erf [r x (k + 2)] + erf [r 2 {k + 2)])(exp [~r x (k + 1)] - exp [-r|(fc + 1)]) 

+ (erf [r x (k)] - erf [r 2 (fc)])(exp [~r x {k + 1)] + exp [-r|(t + 1)])]) 

+ r 3 (fc)[A 4 (fc) - A 2 (fc)](A 4 (fc + l)[(crf [r x (k)} - erf [r 2 (fc)]) 
x (exp [ — 7’j ( k + 1)] — exp [— r| (k + 1)]) + (erf [?’i (fc + 2)] + erf [r 2 (fc -4- 2)]) 

x (exp [ — r x (fc + 1)] + exp [— r 2 (fc + 1)])] ~ A 3 {k + l)[(erf [ri(fc + 2)] + erf [r 2 (Ar + 2)]) 

x (exp [-rfik + 1)] - exp [-r 2 (fc + 1)]) + (erf [r^fc)] - erf [r 2 (fc)J) 

x (exp [-r?(fr + 1)] + exp [-r\{k + 1)])]) 

+ r 3 (k + l)[A 4 (fc + 1) - A 2 (& + l)](exp [~rf(k + 2)] +exp [-r 2 (fc + 2)]) 
x (A 4 (fc)[(erf [ri(fc)] - erf [r 2 (fc)]) + (erf [r x (k + 1)] + erf [r 2 (k + 1)])] 

- A 3 (fc)(erf [r x (fc)] + erf [r 2 (fc)])) (A-8) 
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